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This work investigates linear preceding over non-singular linear channels with additive white 

Gaussian noise, with lattice-type inputs. The aim is to maximize the minimum distance of the received 
• 1-^ 

, lattice points, where the precoder is subject to an energy constraint. It is shown that the optimal precoder 

H ■ 

Cd ■ only produces a finite number of different lattices, namely perfect lattices, at the receiver The well- 

known densest lattice packings are instances of perfect lattices, however it is analytically shown that 
the densest lattices are not always the solution. This is a counter-intuitive result at first sight, since 
previous work in the area showed a tight connection between densest lattices and minimum distance. 
Since there are only finitely many different perfect lattices, they can theoretically be enumerated off-line. 
A new upper bound on the optimal minimum distance is derived, which significantly improves upon a 
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previously reported bound. Based on this bound, we propose an enumeration algorithm that produces a 
finite codebook of optimal precoders. 

I. Introduction 

Linear channel models are very common in communications, and they describe several modem 
communication systems. Examples are multiple antenna systems (MIMO) and OFDM transmis- 
sion. MIMO gained enormous attention with the seminal work in [[Tl, showing the increase 
of capacity with the minimum number of antennas at the communicating terminals. However, 
achieving this rate requires a Gaussian alphabet and Waterfilling (WF), where the latter requires 
perfect knowledge of the singular values and the right unitary matrix of the channel. From an 
application point of view, Gaussian alphabets are not practical, and instead discrete alphabets are 
used in practice. For discrete alphabets, WF does not longer achieve the maximal information 
rateR as demonstrated in [|2l. Instead, [0 derived a technique similar to WF, called Mer- 
cury AVaterfilling (MWF), which is the optimal power loading for discrete signal constellations. 
After the eigenmodes of the channel have been accessed, the power loading corresponds to 
applying a diagonal precoding matrix to the data vector. However, it is possible to achieve 
higher information rates than offered by MWF by applying a non-diagonal precoding matrix, 
as demonstrated in [|3J. A numerical technique that attempts to find the optimal precoder for 
maximizing the information rate for a specific MIMO channel was presented in that work, but 
it was not possible to prove the optimality of the technique. Recently, Hi presented a numerical 
algorithm that converges to the precoder maximizing the information rate. The bottleneck in the 
algorithm is the severe complexity arising from computation of the MMSE matrix at every step. 

Once the discrete data alphabet is large and structured, lattice theory naturally comes into 
focus, since the received signaling points will be organized as a lattice. Lattice theory is a 
subject that has long been studied within information theory, with classical works such as [6J- 
[El. In |l6l, [|7]|, it was shown that there exist codes achieving the capacity over AWGN channels, 
whose codewords are instances of lattice points. The work in [8| categorized the overall gain 
of a code as a coding gain and a shaping gain, both determined by the lattice structure that is 

^The information rate between two discrete sequences y and x is the quantity /(y; x) — H{y) — _ff(y|x), where H{-} is 
the entropy operator for discrete sequences, while capacity is maxp(x) /(y;x), where the maximization is over all probability 
mass functions p(x). 
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used. A different application is [9^1, where the two dimensional lattice structure that minimizes 
the quantization error was found. This formulation gives rise to an optimization problem over 
lattices, and the optimal structure is the well-known hexagonal lattice. An application of lattice 
theory, related to our work, is [10|. Therein, lattice theory was used to construct precoders 
for linear channels that increase the minimum distance of the received signaling constellation. 
Namely, it was suggested that the precoder should organize the received points as the densest 
lattice packing. Based on this assumption, precoders with good bit-error-rate (BER) performance 
under maximum likelihood (ML) decoding were derived. Another, more recent, discovery shows 
the relationship between information rate and minimum distance jSl, BH. Namely, the linear 
precoder that maximizes the information rate in the high SNR regime for discrete alphabets, 
also maximizes the minimum distance of the received signaling points. Hence, the precoder 
maximizing the minimum distance is at the same time minimizing the BER and maximizing 
the information rate at high SNRs, which is a very interesting relationship. Thus, searching 
for precoders which maximize the minimum distance of the received signaling points, is a 
fundamental problem both from a practical and theoretical point of view. 

We follow the problem formulation in [ITOl . and study the lattice structures that the optimal 
minimum distance precoder gives rise to. Some initial progress in this direction was made in 
[fm . which investigated two-dimensional real- valued linear systems with a lattice alphabet at 
the transmitter, and showed that the optimal precoding matrix always produces the hexagonal 
lattice at the receiver. The work in [fT2l extended this to complex-valued lattices, and showed that 
the optimal received lattice, when extended to real-valued space, is always the Schlafli lattice 
D4 (also known as the checker-board lattice). This work extends the two-dimensional results to 
arbitrary dimensions. It will be shown that there are only finitely many optimal structures in an 
arbitrary dimension. 

The paper is organized as follows. In Section HI] we define the system model and the problem 
of interest. The lattice theoretic formulation of our problem is given in Section Unl along with a 
set of necessary mathematical tools. In Section |IVl we show analytically which optimal lattices 
that solve our problem. Bounds are developed in Section |Vl while Section |Vl] presents algorithms 
for finding the optimal solution and also suboptimal precoder constructions, based on the results 
in Sections HV] and |Vl Finally, Section IVIII concludes our work. 
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II. Problem Background 



A. System Model 

We start with some notation used throughout this work. Ej-} denotes the expectation operator. 
7^{- },!{■} denote the real and imaginary parts of a matrix, respectively. M, C and Z denote the 
real- valued, complex- valued, and integer- valued numbers, respectively. Similarly, M*^^^, C*^^^ 
and Z^^^^ denote the spaces of real-valued, complex-valued and integer-valued M x N matrices, 
respectively. S^^^ denotes the space of N x N symmetric matrices and S^q'^ denotes the cone 
of positive definite matrices. Z[i] stands for the set of Gaussian integers. M^, and denote 
the A^-dimensional space of real- valued, complex-valued and integer- valued vectors, respectively. 
Further, Iat is the N x N identity matrix, 0^ the all zero A^-dimensional vector, O^vxiv the N x N 
all-zero matrix, (■)^ matrix transpose, (■)* Hermitian transpose, tr(-) the trace of a matrix and 
||x|| = -\/x*x the Frobenius norm of the vector x. The j:th column of a matrix A is denoted by 
SLj and the element at position in A is denoted by Uij. For a vector x, the j:th element is 
denoted by xj. 

The communication model studied in this work is 



where H G C is the channel matrix and F G C is a precoding matrix satisfying the 
energy constraint tr(FF*) < Pq for some Pq > 0. x = Fa is the transmitted symbol vector, 
whereas a is the precoded x 1 column vector comprising uncorrelated symbols, aj G Z[z], 
with a probability mass function p(a) such that E{aa*} = 1^. In other words, the elements aj 
are crafted from the Gaussian integers and they are uncorrelated, which represents an infinite 
QAM constellation. Finally, n ~ CA/'(Oj\f , Ia/) is complex- valued, circulary symmetric, white 
Gaussian noise and y is the received vector. 

The aim of the paper is to design F such that the minimum distance of the received signaling 
points is maximized. First, we transform ([1} into an equivalent real-valued model. Any complex 
valued MxN matrix Ac is isomorphic to a real-valued 2M x 2N matrix A^, with the isomorphy 
function being 



y = 



HFa + n 



(1) 




(2) 
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For A^-dimensional complex-valued vectors x^, the isomorphy function is 

( ^ 

Xc ^ X,, = I (3) 

where x,. has dimension 2N. Applying these transformations to the matrices and vectors in 
we arrive at a real-valued signaling model 

Yr = HrFrSir + 11,.. (4) 

Thus, without loss of generality, we can work with the model in (HI) rather than ([T). Since a can 
be any Gaussian integer vector of dimension A^, the real-valued vector a,, can be any integer 
vector of dimension 2N. Further, it holds that tr(F^Fr) = 2tr(F*F) < 2Po- The precoding is 
now performed over the real- valued domain as x^ = F^a^, where the actual complex- valued 
symbols x to be transmitted over H in ([1]) are obtained from x^ through the inverse of dSj. 
Note that the transformation in © imposes a skew-symmetric structure on F^, which can be 
relaxed when the precoding is performed in the real-valued domain, i.e., F,, can be any 2N x 2N 
real-valued matrix satisfying the trace constraint. This is not true for H,, since it must represent 
a complex-valued channel. However, since our analysis is applicable to general non-singular 
Hr, we will drop the the skew- symmetric constraint on H,.. Thus, by precoding over the real- 
valued domain, performance gains can be expected because there are more degrees of freedom 
in designing F^ than in designing F. Henceforth, we omit the subscript r and assume that all 
variables in dimensions are real-valued, unless stated otherwise. 

Let e = a — a be the difference vector of two data vectors a 7^ a. The squared minimum 
distance is Dj^j^(HF) = mine^O]v l|HFe|p. Thus, the problem of interest in this work is the 
following optimization problem 

maxD^i„(HF) 

F 

subject to (5) 

tr(F*F) < Pq. 

Let H = USV^ G R^^^ be the singular value decomposition (SVD) of H. We study non- 
singular H, thus we can assume that S has A^ positive diagonal elements. Since U is merely a 
rotation of the received signaling points, it has no effect on the minimum distance and can be 
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discarded from the problem formulation. Also, the matrix V can be absorbed into F without 
changing the transmitted power. Thus, equivalently, we consider the following channel model 

y = SFa + n, (6) 

for which the optimization problem to be studied in this paper becomes: 

maxL)^i„(SF) subject to tr(FF^) < Pq- (V) 

F 

Another formulation of this optimization is possible, by observing that the precoder F solving 
this optimization is the same precoder that minimizes the trace for a fixed value of the objective 
function D^i^(SF). We can write 

L)^.^(SF) = min ||SFef = min e^F^S^Fe = min e^Ge, 

where G = F^S^F. A fixed value D^i^(SF) = d of the objective function implies that e^Ge > 
d, Ve 7^ Oat, where equality is achieved for at least one integer vector e. Since the objective and 
the constraint function in (|7]) are homogeneous of degree 2, we can assume that d = I. Thus, 
an equivalent formulation of the optimization is 

mintr(FF^) 

F 

subject to (8) 

e^Ge > 1 Ve G Z^\{0^}, 

where Z^\{07v} is the set of all A^-dimensional integer vectors except the all-zero vector. Yet 
another equivalent way of expressing ([8]) is to maximize the normalized minimum distance 
o?^i^(S, F) = L)^jjj(S, F)/tr(FF^) over F ^ On-xn- For our purposes, the problem formulation 
in dS]) will turn out to be the most convenient, and will be the one studied in this paper. In the 
next section, we formulate dSJ as a pure lattice problem, and introduce the tools from lattice 
theory needed to analyze it. 

in. Lattice-Theoretic Approach 

This section is split into five parts. Section IIII-AI briefly presents basic lattice theory, while 
Section IIII-BI describes the Ryshkov polytope and Section IIII-CI the Minkowski polytope, both 
of fundamental importance for the understanding of our subsequent analysis. Section IIII-DI 
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formulates ([8]) as a lattice problem, while Section IIII-EI gives an overview of famous lattice 
problems and techniques, applicable to our problem, to solve them. 

We will use some terms from convex geometry in what follows. By an A^-dimensional 
polyhedral cone, we mean the set {AiVi + . . . A^vx : > 0, 1 < j < i^, } for K given A^- 
dimensional points vi, . . . , \k- A polytope in dimensions is the intersection of a finite number 
of A^-dimensional halfspaceo i-c, the set of A^-dimensional points {x : iXi + . . . atXat < 
bj '■ 1 < j < M} for given numbers M, j, hj,\<i<N,\<j< M. A face of a polytope 
is the intersection between the polytope and a supporting hyperplanqj of the polytope. If the 
face is one-dimensional, we call it an edge, or in the case when the polytope is a polyhedral 
cone, an extreme ray. 



A. Lattices 

Let L G R^^^ and let the columns of L be denoted by Ii, . . . , Iat. A lattice Al is the set of 
points 

Al = {Lu : u G Z^}. (9) 

In u is an integer vector and L is called a generator matrix for the lattice Al. The squared 
minimum distance of Al is defined as: 

Dj^jn(L) = min ||L(u — v)|p = min ||Le||^ = min e GlG, 

where u, v and e = u — v are integer vectors and Gl is the Gram matrix for the lattice Al- 
The fundamental volume is VoI(Al) = | det(L)|, i.e., it is the volume spanned by li, . . . , \n- 
Let pj denote a lattice point in Al. A Voronoi region around a lattice point is the set 
Vp^(Al) = {w : ||w — pjll < llpfc — w||, £ Al}. Due to the symmetry of a lattice, it holds 
that Vp^(Al) = Pj + Vo^(Al). The Voronoi region around Oat is denoted V(Al). 

As can be seen from the definition of Al, the column vectors Ii, . . . , I^r form a basis for the 
lattice. There are infinitely many bases for a lattice. Assume that L' is another basis for Al. It 
holds that L' = LZ, where Z is a unimodular matrix, i.e., Z has integer entries and det(Z) = ±1 

^By a half-space we mean either of the two parts into which a hyperplane divides a Euclidean space. 

''a supporting hyperplane of a set 5 is a hyperplane that intersects S, such that S is completely contained in one of the two 
halfspaces determined by the hyperplane. 
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Fig. 1. The hexagonal lattice depicted with a geometrical description of the introduced lattice quantities. 



ifTSl . Hence, the generator matrix L' generates the same lattice as L, i.e., Al = Al/ where = 

T T 

denotes equality between sets. Two Gram matrices Gli = Li Li and Glj = L2 L2 are isometric 

T 

if there exists a unimodular Z and a constant c such that Glj = cZ GljZ. Geometrically, this 
means that Li and L2 are the same lattice up to rotation and scaling of the basis vectors. 
From the definition of the different lattice measures, it follows that 

I^L„(Aqlz) = DlJL) (10) 
where Q is any orthogonal matrix. Similarly, VoI(Aqlz) = VoI(Al). 

B. Ryshkov Polytope 

Let Al C be a lattice with a generator matrix L and with D'^^^{L) > A. This can be 
written as an infinite set of inequalities e^Gi^e > A, where e G /{Ojy} and Gl = L^L. Since 
Gl is a symmetric matrix, its dimension is A^(A^ + l)/2, and the infinite set of inequalities are 
linear over the N(N + l)/2 distinct elements of G. By considering the distinct elements in G 
as a vector {gi^i, . . . , gi^N, 92,2, ■ ■ ■ , 92,n, ■ ■ ■ , 9n,n) in M^^^+i)/^, the infinite set of inequalities 
represent an intersection of infinitely many halfspaces in Next, we define [[T4| 

Definition 1: The Ryshkov polytope Tlx is the set Tlx = {G : e^Ge > A, e G Z^/{Oiv}}- 

It is easily realized that any G G Tlx is positive definite, thus Tlx C S^q^. In the vector 
space the set of positive definite matrices S^q^ corresponds to a cone, where Tlx is 
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contained in the interior of the cone. Tlx is a convex and unbounded set, since if Gi, G2 G Tl\, 
then kiGi + ^262 G TZx for ki, k2 > and ki + k2 > 1. Because any positive definite Gram 
matrix G, hereinafter called a "positive quadratic form" (PQF), corresponds to a lattice, the 
Ryshkov polytope contains all Gram matrices of lattices with minimum distance of at least A. 

Since TZx is the intersection of infinitely many halspaces, it could be the case that TZx has a 
boundary that is "curved" and does not represent a polytope. More formally, there could exist a 
point on the boundary of TZx for which there is only one support plane, which intersects TZx only 
at this point. We say that an intersection of infinitely many halspaces P = fl^^^'Hi, is a locally 
finite polytope, if the intersection of P and an arbitrary polytope is again a polytope. Thus, 
a locally finite polytope P contains no curved boundary. The following theorem [14] justifies 
the name "Ryshkov polytope", and plays a fundamental role for the classification of optimal 
precoders that is developed in this work. 

Theorem 1: For A > 0, the set TZx is a locally finite polytope. 
A vertex in TZx corresponds to a form G that is the unique solution to a set of at least A^(A^+l)/2 
linearly independent equations e^Gej = \, j = 1 . . .K, where K > N{N + l)/2. Note that 

T 

if G is a vertex in TZx, then so is Z GZ where Z is unimodular, and thus there is an infinite, 
but countable, number of vertices in the Ryshkov polytope. This observation also implies that 
the vertices can be partitioned into equivalence classes, where the equivalence relation is an 
isometry between two vertices. 

Lattices corresponding to vertices of TZx are named perfect lattices in the literature [[TSl . and 
the corresponding Gram matrices are perfect forms. The next theorem gives another interesting 
property of the Ryshkov polytope [[TSl . which is important for our work. 

Theorem 2: There are only finitely many non-isometric perfect forms in the Ryshkov polytope. 

Hence, although there are infinitely many perfect forms in the Ryshkov polytope. Theorem |2] 
reveals that out of these, only finitely many are non-isometric and correspond to different lattices. 
The non-isometric perfect lattices have been tabulated for all dimensions up to = 8 [[T4|. In 
two and three dimensions, there is only one unique perfect lattice. In four dimensions, there are 
two, in five there are three, and in dimension = 8 there are 10916. 

Voronoi's algorithm, is commonly used to traverse the vertices [|T5l of the Ryshkov polytope. 
In Section |VI] we present an algorithm for solving our problem, which in essence is a modified 
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version of Voronoi's algoritm. 

C. Minkowski Polytope 

Another characterization of PQFs is via Minkowski reduction. 

Definition 2: The Minkowski reduction region Ai is the set of all G satisfying 

(i) v^Gv > (7j_j, for all v e such that gcd(fi, . . . , vn) = 1. 
{ii) gi,i+i>0, i = l,...,N-l. (11) 

As with the Ryshkov polytope, is a subset of S^^^. A PQF G = L^L G 7W is said to be 
Minkowski reduced and the lattice generator matrix L is called a Minkowski reduced generator 
matrix for Al- It can be shown that any lattice Ab has a generator matrix L that is Minkowski 
reduced, i.e., there exists an L such that B = LZ, where L is Minkowski reduced and Z is a 
unimodular matrix [[T6l . Note that the Minkowski reduced generator matrix is not unique for a 
certain lattice, e.g., if L is Minkowski reduced, then so is — L. However, it can be proved that 
there are only finitely many Minkowski reduced generator matrices for any lattice fT6]. Given a 
generator matrix B, a Minkowski reduced generator matrix L, and the corresponding unimodular 
matrix Z, can both be obtained by applying the Minkowski reduction algorithm on B [16J. 

T 

Let L be a Minkowski reduced generator matrix and Gl = L L the corresponding Minkowski 
reduced PQF. Condition (i) in (fTT)) implies that _D^jjj(L) > gi^i, and since gi^i = it follows 

that L'^;^(L) = Qi i, because at least v = (1 . . . 0)^ achieves equality. Hence, any Minkowski 
reduced generator matrix L contains the shortest vector in the lattice Al as one of its columns. 
Ai is an intersection of infinitely many halfspaces, just as the Ryshkov polytope, but with 
different halfspaces in this case. It is easily seen that Ai corresponds to a cone in the vector 
space M^(^+i)/2, since if Gi, G2 e M , then k^Gi + /tsGa G for fci > and k2 > 0. We 
now define 

Definition 3: Mx = {G : G G M, Qia > A}. 
Hence, the Minkowski reduced PQFs in A^a correspond to all Minkowski reduced generator 
matrices of lattices with a minimum distance of at least A. 

A polyhedral cone is a cone with a finite number of flat faces, and is therefore also a polytope. 
A fundamental result by Minkowski is [1I61I 

Theorem 3: is a polyhedral cone in R^(^+i)/2. 
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Theorem [H is of later importance, since the polyhedral structure of M is crucial for the solvability 
of®. 

It follows from Definition [3] that Aix is the intersection of the hyperplane {G : (711 = A} 
with Ai, and from Theorem |3] we conclude that Aix contains a finite number of vertices. The 
vertices are hereinafter denoted as Minkowski extreme forms, and the corresponding lattices as 
Minkowski extreme lattices. Minkowski extreme forms have been tabulated up to dimension 7, 
while for higher dimensions they are unknown since the computational complexity is tauntalizing. 
Compare this to perfect forms in the Ryshkov polytope, which are known up to dimension 8. 
This is due to the fact that enumerating perfect forms is computationally more tractable than 
enumerating Minkowski extreme forms [14J. Ryshkov managed to show that every perfect form 
is equivalent to a form lying on an extreme ray of the Minkowski reduction region Ai [17]. 
Cohn et al., however, showed that there are extreme rays in the Minkowski reduction region that 
do not contain perfect forms [|T8l . Thus, every vertex (perfect form) of the Ryshkov polytope 
IZi can be reduced to a vertex in A^i, but there are some extreme rays in A^i which contain 
PQFs in IZi that are not vertices of IZi. Thus, since the Minkowski reduction region is different 
from the Ryshkov polytope, defining our optimization problem over it can provide additional 
insights to the properties of the optimal solution. 

D. Lattice-Based Problem Formulation 

We are now ready to reformulate dS) as a pure lattice optimization problem. This will, for 
completeness, be done over the Ryshkov polytope as well as over the Minkowski reduction 
region. We begin with the former. Start by factorizing F as F = S^^UB, where U is an 
orthogonal matrix and B is any matrix such that F satisfies the trace constraint tr(F^F) < Pq. 
From lattice theory, it follows that B can be regarded as a generator matrix for a lattice Ab. 
Inserting the expression F = S^^UB into ([8]), we arrive at 

mintrfB^U^S^^uB) 

U,B ^ ^ 

subject to (12) 
Gb e 7^l, 

T 

where Gb = B B. The subscript in Gb will be left out when no confusion can arise. 
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Let us now instead turn to the second formulation and formulate ([8]) as an optimization over 
the Minkowski poly tope Aii. We keep the factorization F = S^^UB, but we further factorize B 
as B = LZ, where L is a Minkowski reduced basis of the lattice Ab and Z a unimodular matrix. 
This gives that the F in ([8]) can also be factorized as F = S^^ULZ. Furthermore, the constraint 
D^jjj(B) > 1 is now equivalent to Gl = L^L G Aii. Thus, ([8]) can as well be formulated as 

min trfZ^L^U^S^^ULZ) 

U,L,Z 

subject to (13) 
Gl e Ml. 

Also for (fT3l) . the subscript in Gl will sometimes be left out. Note that the objective functions 
in (fT2)) and (fT3l) are exactly the same, since B = LZ, but the optimization procedure is 
different for the two problems. First of all, the optimization domains are different. Secondly, 
the optimization in (fT2l) only involves a minimization over orthogonal (U) and invertible (B) 
matrices, while (fT3l) is a minimization over orthogonal, invertible (L) and unimodular matrices 
(Z). Hence, the methodology for solving (fT2)) differs from the one solving (fT3l) . The formulation 
in (fT3l) also reveals the fact that changing the basis in Ab, i.e., varying Z, only affects the 
transmitted power, which is not evident from the formulation in (fT2)) . Another advantage of 
the formulation in (fT3T ) will be revealed by Theorem |4] in Section Illl-Ci The work in [10| 
considered a problem formulation similar to (flJl . but without using the Minkowski reduction 
domain. Instead, only the objective function was studied for different lattice bases L. Both |8l 
and [[Toll made approximations to the minimum distance problem, the hypothesis was that the 
densest packing lattices in high dimensions should produce large distances. However, no exact 
results were presented. Instead, in [ilOJ it was just proposed that L should be a basis for the 
densest lattice packing, and a heuristic, iterative algorithm was given to find the optimal Z. The 
derived precoders turn out to have good performance, however the question remains whether 
they indeed are optimal minimum distance precoders, and if not, how far away they are from 
the optimum. Thus, (fT3l) was not satisfactory treated in [lOJ. 

Two fundamental questions arise about the problems (fT2] ) and (fT3l) : 1) Is there an explicit 
formula for the optimal solutions to any of the problems? 2) If there is no such formula, what 
is the structure of the solution and can it be found in a simple way for any channel outcome S? 

Before answering these questions, we look at some classical lattice problems and tools for 
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solving lattice optimization problems. Our motivation for surveying these well known problems 
is that the problem studied in this paper is tightly connected, and has in principle the same 
structure in its solution as some of the classical problems. 

E. Classical Lattice Problems 

There are many optimization problems that can be interpreted as optimization over lattices. 
A famous one is finding the densest lattice packing of spheres in an A^-dimensional space, 
corresponding to the following optimization 

min Vol(L) 

^ (14) 
subject toD^;jL)>l, 

i.e., to find, among all lattices with fixed minimum distance, the lattice with the minimal volume. 
The dual of this problem is to find the lattice maximizing the volume of the sphere encompassed 
by its Voronoi region; this is known as maximizing the covering of the lattice. Mathematically, 
it corresponds to the following optimization 

min max ||w|| 

L wGV(L) 

subject to 

(15) 

dLM > 1. 

The general solutions of these problems remain unknown as of today. However, for small enough 
dimensions, solutions are known. In two dimensions, it turns out that the hexagonal lattice solves 
both of these problems; this fact was shown for ([14]) by Lagrange in 1801 fTOl, and for (fTST ) by 
Kershner in 1934 lfT3]| . The packing problem has been solved for < 9 and N = 24, while it is 
unsolved for all other A^. For the covering problem, the solution is known for A^ < 5. Although 
the packing problem is unsolved in general, it is known that the optimal lattice must be a perfect 
lattice which also immediately implies that it is attained at a Minkowski extreme lattice. Hence, 
finding the densest lattice packing in any dimension A^ amounts to traversing the non-isometric 
vertices in TZi, or traversing the vertices in A^i. Although the former is computationally more 
feasible, traversing the non-isometric perfect forms also becomes computationally inefficient for 
higher dimensions. Albeit the computational bottleneck, it is known that (fT4l) and (fT5l) are both 
discrete optimization problems rather than continuous ones. 
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To show that the solution of (fT4l) is achieved by a perfect lattice (or a Minkowski extreme 
lattice), it suffices to show that Vol(L) is a strictly concave function over Sy^^^. Since S^q^ 
contains the polytopes TZi and Aii, this therefore implies that Vol(L) is concave over both TZi 
and Ail. Therefore, the solution to ([14] ) is attained at the vertices of these polytopes, i.e., at the 
perfect lattices (vertices of TZi) and the Minkowski extreme lattices (vertices of A^i). Hence, 
concavity of the objective function is enough to conclude that perfect lattices solve a given lattice 
optimization problem. The concavity of det(G)^/^ over S^^^ was shown by Minkowski [|T6ll . 

Next, we show that the objective functions in (fT2l) and (fT4l) are of different nature. The 
orthogonal matrix U minimizing the objective function in (fT2l) has been found in fTOl, and is 
shown to be equal to the left orthogonal matrix in the SVD decomposition of B. Hence, inserting 
this U into the objective function gives the optimization 

N 

j=i (16) 



subject toGe e TZi, 

where (^^(Gb) is the j:th largest eigenvalue of Gr and Sj is the j:th largest diagonal element 
in S. The optimization in (O can be performed over the Ryshkov polytope, with the objective 
function Vdet(B^S-2B) = ^det(S-2) det(GB). The matrix S can be regarded as a constant 
and does not impact the optimization. It further holds that 



Vdet(S-2) det(Gi 



N 



n^.(GB)A^ 



Hence, (fT4l) minimizes the A^:th root of the product of the eigenvalues of Gb over TZi, while (fT6l) 
minimizes a weighted sum of them. Due to the arithmetic-geometric mean (AM-GM) inequality, 
we have that 



N 



5^c.,(G)/.?>iV 



N 



\ 3 = 1 



which shows that the L solving (fT4l) is only minimizing the lower bound to the objective function 
in (fT6l) . thus not guaranteeing that it is the optimum to (fT6lP Hence, although (fT2l) and (fT4l) have 

'This same reasoning is used in IIIOI in order to propose densest lattices as good candidates for providing a large minimum 
distance. 
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the same optimization domain, (fT2l) posseses a different objective function than (fT4l) . and is thus 
a different lattice optimization problem. 

IV. Optimal Lattice Structure 

This section will prove the concavity of the objective functions in (fT2l) and (fT3l) . respectively. 
We start by proving the concavity of the objective function in (fT3l) over S^^^, for any given 
S and Z matrix. Define /(L, Z) = minu tr(Z^L^U^S^^ULZ), which is the objective function 
in (fT3l) without the minimization over L and Z. Observe also that /(B,lAr) is the objective 
function in (fT2l) without the minimization over B. We now show 

Theorem 4: For a fixed Z, /(L, Z) is concave over 5^0^^ with respect to Gl- 
Proof: Write minu tr(Z^L^U^S-2ULZ) = minu tr(S-2ULZZ^L^U^). Let LZZ^L^ = 
QDQ^ be the eigenvalue decomposition of LZZ^L^, where Q is the orthogonal matrix. Now 
note that 

/(L,Z) = mintr(S"^ULZZ^L^U^) 

= mintr(S-2UQD2Q^U'') 

= mintr(S"^UQ^D2QU^) 

= mintr(S"^UZ^L^LZU^) 
u 

= mintr(S"^UZ^GLZU^). 

Hence 

/(L,Z) = /i(Gl,Z) 

= mintrfS-^UZ^GLZU^). 

Now it follows that for positive semidefinite Gi, G2 and < 7 < 1, 

^(7Gi + (l-7)G2,Z) = min7tr(S-2UGiU^) + (1 -7)tr(S-2UG2U'') 

> 7mintr(S"2UGiU^) + (1 - 7) mintrfS^^^UGaU^) 

= 7/i(Gi) + (l-7)MG2), 
which shows that /i(Gl, Z), and thus also /(L, Z), are concave over S^q^ with respect to Gl- 
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An immediate corollary is 

Corollary 1: /(L, Z) is concave over S^q^ with respect to the Gram matrix Gb = Glz = 
Z^L^LZ. 

Proof: Let L = B, Z = Ijv and apply Theorem |4l ■ 
Taken together, Theorem |4] and Corollary 1 show that the objective functions in (fT2l) and (fT3l) . 
respectively, are both concave over their corresponding domains. This immediately implies that 
the solution to (fT2l) is a perfect lattice, and the solution to (fT3l) is a Minkowski extreme lattice. 
Exactly which perfect lattice/Minkowski extreme lattice that solves ([8]) depends of course on 
the channel outcome S, and an algorithm is given in Section |Vl] that enumerates all possible 
perfect forms solving ([8]) for a specific S. However, since there are finitely many perfect 
lattices/Minkowski extreme lattices in dimensions, we know that there are finitely many 
different lattices solving the problem for all S. This answers our second question posed in Section 
IIII-D[ The optimal Gb in (fT2)) is a vertex of the polytope TZi, and the optimal Gl solving (fT3l) 
corresponds to a vertex in A^i. Hence, the solution to ([8]) does not depend continuously on 
S, instead it changes in a discrete fashion when S is varied continuously. Relating to our first 
question in Section IIII-Dl this result implies that an explicit formula for the solution of ([8]) 
seems out of reach, since such a formula does not exist for (fT4l) whose set of possible solutions 
is a subset of the set of possible solutions to ([8]). Altogether, a previously unknowned result 
is revealed: There are finitely many lattices that can solve the minimum distance optimization 
problem in ([8]), and they can be enumerated offline. 

Theorem |4] also reveals that for any given Z matrix, the optimal solution to (fT3l) occurs at 
a Minkowski extreme lattice. Thus, given any Z, the optimal L that builds up F = S^^ULZ 
in dH]) is a Minkowski extreme lattice. This fact will be used in Section |Vl] to develop a good 
suboptimal precoder construction. Hence, the problem formulation in (fT3l) provides additional 
information about the behavior of (|8]), not present in (fTZl i: This is the main reason for introducing 

We can already at this stage deduce several interesting conclusions from the result in Theorem 
m For up to three dimensions, there is only one non-isometric perfect lattice in each dimension: 
For A^ = 2 it is the hexagonal lattice and in A^ = 3 it is the face-centered cubic lattice. Since 
these are the only non-isometric perfect lattices in these dimensions, they also solve (fT4l) . and 
thus the proposition in [[TOl to use densest lattice packings in (fT3l) is optimal for these dimensions. 
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However, when = 4 there are two non-isometric lattices: The checkerboard lattice D4 and 
the root lattice A4 [fT3l. It will be demonstrated in Section IVI-CI that both of these lattices occur 
as solutions to (fT2l) for different S, so the constructions in [[TOl are suboptimal for = 4. 
Another interesting consequence of Theorem |4] is that the main result in (fTT], which shows that 
the hexagonal lattice is optimal in two dimensions, now follows immediately from Theorem |4] 
Note that it is now an easy task to construct the optimal F in once the Gb solving 
([I2I) is known. Let Gopt = B^p^Bopt denote the optimal form and Gopt = QoptDoptQopt be its 
eigenvalue decomposition. Since the optimal U in F = S~^UB is equal to the left orthogonal 
matrix in the SVD decomposition of B, it follows that the optimal F can be constructed as 

Fopt = S-iv^U;;pt. (17) 

To summarize, we now have the following knowledge at hand about the solution to our original 
problem in ([8]). We have shown that ^ is equivalent to both (fT2l) and (fT3l) . Theorem |4] then 
shows that the L matrix solving (fT3l) . for any invertible S, gives rise to a Gram matrix G = L^L 
that corresponds to a vertex in the polytope Aii. Since Aii has a finite number of vertices for 
any dimension A^, and is independent of the matrix S, it holds that there are finitely many L 
matrices (up to rotation) that are candidates to solving (flSl) for any given S. Once the optimal 
Lopt is known (up to rotation), it remains to find the optimal unimodular matrix Zopt in (fT3l) 
and then to construct the optimal precoder Fopt from (flTI) . where Bopt = Lopt Zopt- To find the 
optimal Z, we need to perform a search over unimodular matrices, which can be simplified if 
good bounds on the optimum solution to (fT3l) are known. These bounds will be developed in 
the next section. 

When it comes to the equivalent problem formulation in (fT2]) . Corollary [T| shows that the B 
matrix solving (fT2]) for any given invertible S, is such that it produces a Gram matrix G = B^B 
that is one of the vertices in the polytope TZi. Given the optimal G in T^i, the optimal precoder 
is obtained through (fTTT ). The TZi polytope contains infinitely many vertices, and it is known 
that each vertex is isometric to some vertex in A^i. Since (fT2)) is connected to (fT3l) through 
the factorization B = LZ, it holds that if Gb = B^B = Z^L^LZ is a vertex in TZi, then 

T 

Gl = L L is a vertex in A^i. Hence, traversing the different vertices in TZi is equivalent to a 
joint enumeration of some of the vertices in A^i and different unimodular matrices Z. However, 
it is instead possible to directly enumerate perfect forms by formulating an algorithm working 
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over TZi. Again, bounds are needed in order to restrict the amount of vertices to traverse, and 
they will be presented in the next section. 

V. Bounds On the Optimal Solution 

We start by deriving lower and upper bounds to tr(Z^L^U^S^^ULZ), which is the objective 
function in (fT3l) . The upper bound presented here improves significantly upon the upper bound 
presented in [10]. From these bounds, we are able to derive further bounds that aid in restricting 
the search space for the algorithms to be introduced in Section |VIl 

Theorem 5: The following lower bound holds for the optimal solution to (fT3l) 



min tr(z'^L''u^S"2ULZ) > ^<ydet(L)/ det(S). (18) 

XJ.L, Z 



Proof: Dropping the integer- valued constraint on Z, while keeping the determinant constraint 
det(Z) = ±1, we apply the method of Lagrange multipliers to find first order optimality 
conditions. Let M = L^U^S^^UL. The optimal Zo must satisfy 
9tr(Z^L^U^S-2ULZo) det(Zo) - 1 



where 7 G M. Taking determinants on both sides, and making use of det(Zo) = ±1, we get 
7 = 2 ^ det(M). Inserting this 7 into (fT9l) and multiplying both sides of the equation with Z^, 
we arrive at Z^MZq = ^det(M)Ijv. Hence, for this Zo, we get 

tr(ZoL''u^S~2ULZo) = tr(I^) ^det{M) 

= N^det{M). 

Since this expression is independent of U, it is a lower bound to the objective function for a 
fixed L. Expressing det(M) = det(L^S"^L) = det^(L)/ det^(S), we arrive at the lower bound 
in ([H. ■ 
This bound was also reported in [|T0] . but derived in a different way, by using the AM-GM 
inequality and Hadamard's inequality. The approach presented here shows that this lower bound 
corresponds to the optimal real-valued unimodular Z. 

Next, we derive an upper bound on minu,L,z ti'(Z^L^U^S^^ULZ) 

Theorem 6: If D^^^(L) = 1, then 

min tr(Z^L^U^S"^ULZ) < A^ ^/yiTdetTS). (20) 

U,L,Z V / \ / 
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Proof: Let B = ULZ = QR denote the QR-decomposition of the received lattice. It holds 
that D'^^J^B) > mini<j<jv Ir^^ip lEOj. In lEH, an orthogonal precoder matrix Fgmd (geometric 
mean precoder) and an orthogonal receiver matrix Wg^d were constructed, such that in the QR 
decomposition of Bgmd = WgmdSFgmd, all diagonal elements of R equal ^det(S'). Hence, in 
essence, the precoder Fg^d together with the rotation Wg^d at the receiver, produces a lattice with 
maximal lower bound on D'^-^^. This value is equal to the geometric mean of its singular values, 
thereby its name the "geometric mean precoder". It is clear that D^j^(Bgmd) > v^det(S'), and 
since Fg^d is orthogonal, tr(Fg^dFg„,d) = N. Hence dl^^^{S, Fg^d) = D'^-^{Bg^d)/tT{Fl^^Ff^^d) > 
^^^det{S)/N. Now it follows that for any precoder F with higher (i^in(S, F) than (i^in(S, Fg^d), 
<ij,(S,F) > (det(S))2/^/A^. Hence, this gives an upper bound on tr(FF^), tr(FF^) < 
^/^l/det(S). Writing F = S^^ULZ, we get the upper bound in (EQ]). ■ 
Combining Theorem |5] and |6l we have the following bounds 

iV(det(L)/det(S))2/^ < tr(Z^L''u^S-2ULZ) < iV(l/ det(S))2/^. (21) 

(f2T)) translates into 

(det(S)/det(L))2/^/iV > rfL„(S,F) > (det(S))^/^/iV. (22) 

Note that the L in and ([221) is such that -Dmin(L) = 1- Also, these bounds hold for any 
precoder F = S^^ULZ that improves upon Fgmd- It is readily seen that the ratio between the 
upper bound and the lower bound in (|2TI) is (1/ det(L))^/^. This can be compared to the bounds 
in [[Toll , where the ratio between the upper bound and lower bound contains the exponential factor 
2Af/2 jj^us, the improvement in the upper bound is significant. Hence, for a fixed dimension, the 
optimum ratio (i^in(SFopt)/(det(S))^/^ is always smaller than (l/det(L))^/^, independently 
of the channel S. For example, when N = 2, the optimal lattice is the hexagonal lattice Lhex 
and when D^^^(Lhcx) = 1, det(Lhcx) = 1/2- The ratio between the upper bound and lower 
bound in (|22)) is then 2, hence the optimal (i^i„(SFopt) in two dimensions is at most twice the 
lower bound dl^-^^^SFg^d) = det(S)/2. For = 3, the optimal lattice is the face-centered cubic 
lattice L^a that has a volume of det(L^3) = 1/2 when D^^^(Las) = 1- In this case, the ratio 
between the bounds in (f22)) is 2^/^ ^ 1.59; hence, the performance of the optimal precoder is at 
most 59% better than for Fgmd- It is worthwile to observe that the ratio between the bounds, for 
optimal packing lattices L, equals Hermite's constant 5n = maxL -D^in(L)/Vol(L)^/^. Hermite's 
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constant is the ratio between the constraint function and the objective function in (fT4l) . and is 
therefore an optimization problem equivalent to ([HJ. The following upper and lower bounds are 
known for 5m IfBl. Il22ll 

+ ^ + Cna <5n< — (1 + Cn,2) 



2716 ' 27re ' '"'^ — ^ _ ^ ' ^^-^ ' 

where cat i and C7v,2 are constants depending on the dimension A^. From this it follows that Sn 
grows linearly with the dimension A^, and thus the ratio of our bounds grows at most linearly 
with A^. Note however that Fgmd operates above the lower bound in (|22]) . and due to its good 
SER performance as reported in [|2T| . it can serve as a basis for developing a suboptimal Z 
matrix to (fT3l) . This will be presented in Section |Vll 

As discussed in the first paragraph of Section |IVl a closed form solution to ^ seems out of 
reach. We are thus interested in an algorithm that can find the optimal Z in (fT3l) . or an algorithm 
to find the optimal Gb in (fT2l) . In order to do so, it is desirable to first have some bounds on 
the Z matrix or some quantity depending on it, in order to restrict the search space. From this 
perspective, we develop an upper bound on tr(Z^L^LZ). 

Theorem 7: With Z}^jj^(L) = 1, the following upper bound holds 

mintr(Z^L'^LZ) < N [ , \ . (23) 

Proof: Let Gl = Z^L^LZ. Inserting the optimal U into (|2TI) . we arrive at the upper bound 

tr(^](GL)S-') < iV(l/det(S))2/^, (24) 

where ^^(Gl) is the diagonal matrix containing the eigenvalues (X'j(Gl) of Gl- Since the 
eigenvalues uoj{Q<\) are sorted in opposite order to sj^, and < ■ ■ ■ < s]^^, we have the 
inequality 

which gives us the upper bound in (|23]) ■ 
Geometrically, the inequality in (|23l) implies that the lattice vectors of the optimal lattice LoptZopt 
must have a bounded length. Using the trace inequality [|23l 

T T T T 

|^Ar(GLopt)tr(ZoptZopt) < tr(ZQp^LQptLoptZopt) < a;i(GLopt)tr(ZoptZopt), 
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we also have the following upper bound for Z, 



'opt 



(Zoptz;j < ^^^^^^^^^ ^^/^ j 



2 



(25) 



In terms of Gb = Z L LZ = B B, the upper bound in (|23T ) is 



) 



2 



tr(GB) < N 



ydit(s) 



(26) 



Let ub(S) denote the upper bound in (|26] ). Hence, the optimal Gb in the Ryshkov polytope that 
solves (fT2l) is one of the vertices of the finite, bounded polytope TZi fl {Gb : tr(GB) < ub(S)}. 



In this section, we present algorithmical approaches to solve (fT2l) and (fT3l) . along with a novel 
suboptimal precoder construction based on Theorem |4] Beside this, by applying the knowledge 
that perfect forms solve (fT2] ). we provide a numerical example where the densest lattice packing 
is not a solution to (fT2l) . This shows that the solution to ^ is somewhat counter-intuitive: The 
optimal packing of points at the receiver, does not always minimize the total energy of the lattice 
points at the transmitter. 

A. Methods to solve (O and (fT3l) 

Although (fT2l) can be solved by enumerating all vertices in the polytope T^i fl {G : tr(G) < 
ub(S)}, the methodology for solving (fT3l) will provide another interesting observation. Addi- 
tionally, the problem formulation in (fUl ) gives novel insight into efficient suboptimal precoder 
construction, not present in the formulation in (fTZl i. We divide this section into two parts: First 
we discuss a method to solve (fT3l) . then we discuss the solution to (fT2l) . 

1) Finding the solution to (fT3l) ." To find the solution to (fT3l) . one needs to tabulate Minkowski 
extreme lattices in dimensions. Unfortunately, it turns out to be more complex to enumerate 
Minkowski extreme lattices than perfect forms [14] . Note, however, that only those Minkowski 
extreme lattices that correspond to perfect forms have to be known. Namely, once all non- 
isometric perfect forms have been tabulated in A^ dimensions, the fact that each perfect form 
(vertex) in IZi is equivalent to a Minkowski extreme form (vertex) in M.i, implies that the 
Minkowski extreme lattices solving (fT3l) are the ones corresponding to the non-isometric perfect 



VI. Numerical Methods 
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forms. Therefore, it is not necessary to know all the Minkowski extreme lattices in dimensions 
in order to solve (fTJt . only those corresponding to perfect forms are needed. However, when 
constructing a good suboptimal solution to (fT3T ), presented in Section |VI-B[ it is necessary to 
know all the Minkowski extreme lattices to obtain the best suboptimal construction. The smallest 
eigenvalue in (l25l) is non-zero for all the Minkowski extreme lattices L (candidates for 

the optimum), and thus the bound in (l25l) is well-defined. A geometrical interpretation is that 
this inequality bounds the squared lengths sum of the basis vectors in the integer lattice Z^, 
where the basis vectors are now the rows of Z. Thus, finding the optimal Z can be regarded as 
searching for basis vectors inside a sphere of a certain radius. If one has a priori knowledge about 
the maximum ratio si/ det(S), an off-line, one-shot algorithm can be formulated that searches 
for unimodular Z inside the largest sphere, corresponding to the Minkowski extreme lattice L 
with smallest u;iv(GL) and the channel S with largest upper bound in (|25T ). This sphere certainly 
includes the optimal Zopt corresponding to the optimal Minkowski extreme lattice Lopt for any 
channel that can occur. A large codebook of matrices LZ can then be constructed off-line, by 
multiplying each encountered Z in the sphere with the different Minkowski extreme lattices and 
storing the resulting matrices into the codebook. To then find the optimal precoder F = S^ULZ 
online for a certain S, one simply goes through every element LZ in the codebook, and constructs 
F by using the optimal U. 

The outlined method to solve (fT3]) includes the following steps: 1) Find all Minkowski 
extreme lattices that correspond to non-isometric perfect forms. This is accomplished by ap- 
plying Voronoi's algorithm to enumerate non-isometric perfect forms [24J, and then applying 
the Minkowski reduction algorithm to the obtained perfect lattices. 2) Enumerate all unimodular 
Z satisfying the bounds in (|25l) . There are specialized algorithms for this task [25] . 

2) Finding the solution to (fT2l) .- As described in Section |Vl the optimal G is one of the 
vertices in the polytope TZi fl {G : tr(G) < ub(S)}. Hence, one method to find the optimum 
is to directly enumerate all the perfect forms inside the polytope. A finite codebook can be 
constructed off-line if a priori knowledge of the upper bound in (l26l) is available. A method 
to enumerate perfect forms is via Voronoi's algorithm [24]. It enumerates perfect forms and 
stops when all non-isometric forms have been found. As mentioned, for today's computers, it is 
only usable up to 8 dimensions due to the large number of edges in the Ryshkov polytope in 
higher dimensions. We need to slightly modify the classical Voronoi's algorithm, by changing its 
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Stopping condition. Since we are interested in forms that are isometric, our stopping condition 
is based on the upper bound in (|26] ). Further, we note that transforming Gb to n^GBH, where 
n is a generalized permutation matrix 11 such that each non-zero element in 11 is either 1 or 
—1, does not change the value of the objective function in (fT6l ). since the eigenvalues of Gb 
are still the same. Moreover, this does not change the constraint region in (fT6] ) since it merely 
permutes the error vectors. For our algorithm, this observation implies that G matrices such 
that G = n Gn, for some already visited perfect form G, do not have to be traversed by 
the algorithm; thus, the search space can be reduced. The algorithm needs an initial perfect 
form as starting position, and a good starting point is the root lattice An [14J. The following 
notation is used in the algorithm. Min(G) denotes the set of minimum vectors of G, i.e., the 
set {x : x^Gx =1} and G[x] = x^Gx. The algorithm is summarized by the pseudo-code in 
Table U The only difference between Algorithm 1 and the Voronoi algorithm presented in [24] 
is the stopping condition and search space reduction through generalized permutation matrices. 
Voronoi's algorithm stops as soon as all neighbouring perfect forms of a certain perfect form 
are isometric to some other perfect form already encountered. Algorithm |I] stops as soon as all 
perfect forms (up to a generalized permutation matrix) satisfying the upper bound in (|26] ) have 
been enumerated. 

The Fincke-Pohst algorithm is used to compute Min(G) [26J. The toughest part of the 
algorithm is to compute the extreme rays of a polytope specified by linear inequalities. This is 
the bottleneck of enumerating non-isometric perfect forms with Voronoi's algorithm and thereby 
solving the lattice packing problem in high dimensions. There exist methods that do this in 
0{Nvd) time, where is the dimension, d the number of non-redundant inequalities describing 
the polytope and v is the number of vertices in the polytope [|27l . 

For step 2, determining the neighbouring perfect forms can be done by the algorithm in 
[[T4| [Algorithm 2, Chapter 3], which computes the a needed in step 2. In the other steps, we 
use boolean variables bc^ to denote whether a vertex has been visited or not. 

B. Suboptimal Precoder Construction 

Finding the optimal solution is a computationally demanding task for today's computers, and 
suboptimal solutions are of interest. We base our suboptimal construction on Fgmd from Section 

m 
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Algorithm 1 

Input: A starting perfect form G^, e.g., the root lat- 
tice An. 

Output: The list of G matrices corresponding to the 

VERTICES in the POLYTOPE H {G : G < ub(S)} 

Let G = Gs and define the boolean variable ba, = 1. Save the 
pair (G, ba) in a sel Q = {(G, ba)}- 

1) Compute Min(G) and the extreme rays (edges) 
Ti, . . . , Tfe of the polyhedal cone 

{G' e S^''^ : G'[x] > Vx G Min(G)}. 

2) Determine neighbouring perfect forms Gi as G, = G + 
aTi, i^l...k. 

3) Let Gjj , . . . Gj„ be those neighbouring forms satis- 
fying the upper bound in i26\ and that cannot be 
expressed as 11 GII for some (G,6q) in Q. De- 
fine &Gjj = 0, I — 1 . . .m and let t/ = Q U 
{(G,,,6g,J,...,(G„„,6g^.^)}. 

4) Find a pair (Gj , ) in Q such that baj ~ 0. If such 
a pair exists, change the value of be. to ba^ — 1, let 
G = Gj and go to step 1. Otherwise, stop and return Q. 



TABLE I 

An algorithm that traverses all perfect forms satisfying the bounds in l |26l l. 



It can be numerically verified that for Fgmd, the received lattice is not a Minkowski extreme 
lattice, and is thereby not optimal. Hence, the performance of Fgmd can be improved by ap- 
plying the result of Theorem |4] Let Bgmd = SFgmd be the received lattice at the receiver, 
where Fgmd is scaled so that D^^^(Bg^d) = 1. Now perform a Minkowski reduction on Bgmd 
by using the Minkowski reduction algorithm [[T6l . so that we can factor the basis matrix as 
Bgmd = UgmdLgmdZgmd for somc rotation Ugmd, Minkowski reduced lattice basis Lg^d with 
D^in(Lgmd) = 1, and unimodular Zgmd- It then follows that Fg^d = S^^UgmdLgmdZgmd- 
Define F^.i = S^^UjLm jZgmd, « = 1 ... -ft', to be the K different precoders where j is 
the i:th Minkowski extreme lattice with D'^^^(Lm,i) = 1, and Uj is the left orthogonal matrix 
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of Lm jZgmd- Applying Theorem |4l we know that tr(F^ ^F^j) < tr(Fgj^jFgind) for some 
1 < j < K- Thus, Frnj is a precoder performing better than the geometric mean precoder. 

Hence, by performing a Minkowski reduction and using the resulting unimodular matrix 
together with one of the Minkowski extreme lattices, it is possible to improve upon the geometric 
mean precoder and reach closer to the lower bound given in (^Th . However, performing a 
Minkowski reduction includes finding the shortest basis vector in the lattice, which is an NP- 
hard problem [|28l . Nevertheless, it turns out to be easily doable with a standard workstation at 
least for < 15. Another method that can be used for this purpose is the iterative algorithm 
presented in [[TOl . 



C. Packing lattices are not always a solution to (|T2l) 

Note that a perfect form G that is a candidate for solving (fT2] ) must satisfy the upper bound 
(|23]). Since e^Ge > 1, Ve G Z^/{Ojv}, it holds that Qi^i > 1, i = 1 . . . A^, and thus tr(G) > A^. 
Let \i denote the z:th shortest vector in the lattice L with Gram matrix G. By definition, the 
minimum distance is Ai = 1. Now assume a channel S such that the upper bound in (|23] ) is 
smaller than (A^ — l)Ai + A2. If Z^GZ is another perfect form isometric to G that is also 
a candidate for solving (fT2]) . then tr(Z^GZ) < (A^ — l)Ai + A2. However, this inequality 
significantly limits the number of possible Z matrices and thus the number of perfect forms 
isometric to G. Namely, since tr(Z GZ) = X]j=i ^J^j ^'^^ each term Ai < Z^ GZj < A2, 
it follows that each Zj must correspond to a minimum vector of G, i.e., Zj belongs to the set 
Min(G). 

Let us apply this idea to 4-dimensional lattices. In 4 dimensions, there are only two non- 



isometric perfect forms, the 1)4 and Aa, lattice. A Gram matrix for D, 

/ 1 0.5 \ 
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1 



(27) 



^Gram matrices for non-isometric perfect forms can be found at 1291 . 
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and for A4, 
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-0.5 








-0.5 
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(28) 



Hence, any perfect form in 4 dimensions can be expressed as either Z G/j^Z or Z Ga^'Z for 
some unimodular Z. Further, it holds that Ai = 1 and A2 = 2 for both G^^ and G^^. Now let 
S be 



/ 1 \ 

0.95 
0.94 



(29) 



\ 0.93 / 

The upper bound in (|23T ) is 4.83 for this S. Hence, if a perfect form Z^G^^Z isometric to G^^ 
solves (fT2l) . then the columns of Z must be taken from Mm{Gr)_^). Similarly, if a perfect form 
isometric to solves (fT2l) . then the columns of the corresponding unimodular Z are taken from 
Min(G/i4). It is an easy task to find Mm{GD4) Min(GA4), by applying the Fincke-Pohst 
algorithm, and also to find all unimodular matrices whose columns consist of these minimum 
vectors. Going through each perfect form obtained from these unimodular matrices, and plugging 
in the optimal precoder (fTTT ) into (fT2]) . the result is that the perfect form isometric to Gyi^ gives 
the smallest value of the objective function in (fT2l) . Repeating the same argument for the channel 

/ 1 \ 

0.99 

0.94 

\ 0.93 / 

one concludes that the perfect form isometric to G^^ solves (fT2)) . Hence, this shows that both 
A4 and D4 occur as optimal lattice structures at the receiver; which one it is, depends on the 
channel S. 



(30) 



VII. Conclusions 

This work studies precoding over non-singular linear channels with full CSI through a lattice- 
theoretic approach. The classical complex-valued linear channel is first transformed to a more 
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general real-valued model which enables performance improvements over the classical complex- 
valued model. Then, the main problem studied in the work is to find lattices that maximize the 
minimum distance between the received lattice points, under an average energy constraint at the 
transmitter. The optimal lattice is analytically shown to be a perfect lattice, as defined by Ryshkov, 
for any given non-singular channel. Bounds on the optimal performance are developed, tighter 
than previously reported, which enable construction of algorithms that produce a finite codebook 
of matrices, from which the optimal precoder can be derived. Furthermore, a suboptimal precoder 
construction is presented together with bounds on its performance, which is analytically shown 
to improve upon a previous presented precoding scheme in the literature, by utilizing the new 
results in this work. In addition to this, we demonstrate with an example that optimal packing 
lattices are not always optimal for maximizing minimum distance, which is a counter-intuitive 
result at first sight. An immediate practical application of the derived results is precoding over 
large alphabets. 
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